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We study in this paper the edge polarizations and their consequences for a biased Bernal stacked 
bilayer graphene nanoribbon with zigzag termination. The magnetic states are classified according 
to the interlayer and intralayer couplings between the edge polarizations, and the magnetic phase 
diagram of doping versus bias voltage is given. Coplanar magnetic phase is found and the variation 
of its magnetic structure with the bias voltage is investigated. For all the magnetic states, we also 
discuss the possibility of the half-metallicity, and for a ribbon with perfect zigzag edges we predict 
seven kinds of the half-metallic states, which are characterized by their distinct magnetic structures 
and quantized electrical conductances along the ribbon. 
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I. INTRODUCTION 

Bilayer graphene has attracted considerable attention 
due to its unusual electronic properties and prospective 
applications in graphene-based spintronics 1 -^. While the 
low-energy excitations in monolayer graphene are mass- 
less Dirac fcrmions with linear dispersion and carry Berry 
phase 7T— , the chiral quasiparticles in Bernal stacked bi- 
layer graphene show finite mass with parabolic disper- 
sion and carry a 2w Berry phased. For bilayer graphene 
at charge neutral point, different to the gapless feature 
of single-particle excitations, theories 5 - - — recently predict 
some broken symmetry states due to strong electron- 
electron interactions, which are characterized by small 
finite gaps, supported by some recent experiments 1 ^ - — . 

Another strategy for achieving an interaction-induced 
gap for bilayer graphene is to construct a nanoribbon, 
especially with zigzag terminatio n 15 ' 16 , where the gap 
opening is due to the interedge interactions between the 
spin polarizations of the localized edge states. The edge 
magnetism of zigzag-terminated graphene nanoribbon is 
of great interest for many scientists even though strong 
experimental evidence for its existence is still absen t 17 ' 18 . 
According to the Stoner criterion, at or near half-filing, 
the ground state has a ferromagnetic (FM) instability 
due to electron-electron interactions and will then lead 
to an antiferromagnetic (AF) structure where each edge 
is ferromagnetically polarized but the edges are coupled 
with each other antiferromagneticall y 19 ' 20 . In the low 
carrier doping regime, it was predicted that besides the 
collincar magnetic states, the monolayer zigzag graphene 
nanoribbon(MZGNR) can have a noncollinear(NC) mag- 
netic ground stated—. 

Bilayer zigzag graphene nanoribbons(BZGNRs) under 
a bias voltage provide an additional way to tune the 
electronic properties of graphene-based spintronics^ 4 - - — , 
since charge transfer between the layers affects dra- 
matically the bulk and edge states. Thus the bilayer- 
graphene-ribbon-based electronic devices with the mag- 



netic and transport properties tuned by an electric filed 
can be expected 2 ^—. To the best of our knowledge, pre- 
vious works on the edge magnetism in BZGNR mainly 
focused on the undoped half-filling case and only the 
collincar spin polarizations were taken into account^! - — . 
Therefore the possibility of stable NC magnetic states 
in a doped BZGNR becomes a natural question which 
deserves further investigations. 

On the other hand, under the application of a 
transverse field across a nanoribbon, transition be- 
tween the semiconducting and half-metallic(HM) states 
can be achieved by tuning the field strength in both 
MZGNRi&ii and BZGNR 4 - 1 , For the MZGNR, in the 
absence of the inversion symmetry, which may be caused 
by the substrate, it is found that the HM states can be 
stabilized at both the charge neutrality poin t 45 ' 46 and fi- 
nite doping^ 5 -. In analogy with MZGNR, half-metallicity 
can also be expected in a biased BZGNR due to the ex- 
istence of edge states and their edge polarizations 4 ^. 

In this paper, based on the tight-binding model for 
a biased BZGNR and taking the Coulomb interactions 
into account, we theoretically explore and classify all the 
possible magnetic states due to the couplings between 
edge polarizations and then investigate the possibility of 
the NC mangnetic states. We also search the possibil- 
ity of the half-metallicity for all the magnetic phases. It 
is found that the NC states do exist and can be stabi- 
lized as metallic states in a lightly-doping and small-bias 
regime. Several kinds of HM states are found, which are 
distinguished by their unique magnetic structures and 
quantized electric conductance along the ribbon. 

The paper is organized as follows. In section 2 we give 
the model Hamiltonian and its mean- field treatment. In 
section 3 we classify and discuss the magnetic structures 
of the ground states and then give the phase diagram. In 
section 4 we investigate properties of the NC spin canted 
states and find out all the possible HM states. In Section 
5, we summarize our results. 



II. MODEL 

Different BZGNRs are distinguished by the inter- 
layer stacking arrangements and the edge alignments. 
Here we study the Bernal stacked BZGNR with the /3 
alignmen t 37 ' 38 . It consists of two coupled graphene rib- 
bon sheets, where the sublattice A sites of the upper layer 
are located exactly on the top of sublattice B sites of the 
lower layer. We assume that the system is subjected to 
an external electric field between the layers. The model 
Hamiltonian can be given as follows, 



H = -t J2 (<4,,:<x c tj<x + h.c. 
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where t is the intralayer nearest-neighbor hopping inte- 
gral, and c Ti i a is the electron annihilation operator at 
site i with a =\, \. and r — ± denoting the upper and 
lower layers. In the second term, the interlayer coupling 
t± is the hopping integral between the B sites on the 
lower layer and their corresponding top A sites on the 
top layer. U is the on-site Coulomb repulsion energy, and 
2V the electrostatic bias voltage. Despite the fact that 
V term breaks the inversion symmetry, the Hamiltonian 
at half-filling is still kept invariant under a combination 
operation of PI, where P is the particle-hole transforma- 
tion c Ti icr— >ric* Tia , with 77 = 1(— 1) if i belongs to A(B) 
sublattice of the upper or lower layers) and / the space 
inversion, with the inversion center chosen as the center 
of the bilayer zigzag ribbon. 

Now that the possibility of the noncollinear spin po- 
larization can not be excluded, the Hubbard term is so 
decoupled that the mean-field Hamiltonian can be writ- 
ten as, 
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where the electron spin polarization and charge density 
are given by m Ti = (l/2)(cl l crc T:i }, and n Tti = (c|-c Tii ), 
respectively, with er = (cr x ,<7y,a z ) the spin Pauli matri- 
ces and c' ri = ( C T i+ , C T i\) ■ The last constant term is 
E Q = -U'Y, iT (n^j4 - ml ti ). We choose t w 2.8eF, 
U/t = 1.0 and t±/t = 0.1, which are the established 
values for the coupling parameters^ 8 - - — . For this inter- 
mediate U , the decoupling process introduced is believed 
to be reliable and effective. 




FIG. 1: (Color online) Band structure for a BZGNR with 
width W in the absence of Coulomb interactions (U = 0): 
(a)V/i = 0, W = 200; (b)V/t = 0.1, W = 200; (c)V/t = 0, 
W = 10; (d)V/t = 0.1, W = 10. 2V" is the electrostatic bias 
voltage. The width of a BZGNR W is defined as the number 
of zigzag chains per layer and energy is measured in unit of 
t = 2.8eF. 



We consider a strip bilayer graphene sample with 
zigzag termination which is periodic along the vertical 
direction but open along transverse direction. Equation 
(2) can be diagonalized to be H = J2 n Enlnln + E n 
by introducing quasiparticle annihilation operators j n 
via c T 



= E r 



T 7„, with u™ ia the expansion coef- 



ficients or amplitudes of the eigenstate corresponding to 
eigenvalue E n . At zero temperature, the self-consistent 
parameters are given by n T ^ = EniEtrl^wP an d 

m ra = YsnUi 1 / 2 ) T, a p K^^K^. Here, to get dop- 
ing dependent properties, we fix the total electron num- 
ber 7V e rather than the chemical potential of the system. 
Starting from different initial random spin and density 
configurations, several self-consistent magnetic solutions 
may be found by the standard iteration method. By com- 
paring their total energies (J2n=i E n ) + Eq among these 
solutions, the ground state can be determined by finding 
out the lowest-energy one. In the following discussion, 
the width of a BZGNR W is defined as the number of 
zigzag chains per layer, while the effective doping value is 
defined as Sn = (1 — n) x 2W with n the average electron 
number per site. 



III. THE COLLINEAR MAGNETIC 
STRUCTURES AND THE PHASE DIAGRAM 

In the absence of the Coulomb interactions between 
electrons, an unbiased BZGNR has four partly flat bands 
at E = which correspond to the four edge states with 
two per edge, while for a biased BZGNR, the four de- 
generate flat bands are split into two families, each of 



which consists of one fiat edge band and one nearly flat 
edge bancU^. The band structures in these two situa- 
tions are demonstrated in figure [TJ Therefore, at or near 
half-filling, according to Stoner criteria for magnetic in- 
stability, when the Coulomb interactions are taken into 
account, no matter how small they are, edge polariza- 
tions are expected to occu r 19 ' 51 ! 52 . For a BZGNR with a 
width W = fO, by an extensive numerical investigation 
on model Hamiltonian (2), we find many magnetic states 
with edge spin polarizations. These magnetic states con- 
sist of fourteen types of collinear configurations(figureE|), 
as well as the NC canted states including both the sym- 
metric and asymmetric ones. In this section we focus on 
the collinear states and leave the discussion on the NC 
ones to next section. Similar to MZGNR, one can classify 
the fourteen collinear states into five categories, accord- 
ing to the intralayer and interlayer couplings between the 
edge polarizations. 

In figure [3l the magnetic phase diagram for a bi- 
ased BZGNR with width W = 10 is given. In the 
light-doping and small-bias regime, both edges are spin 
polarized. When both the edge spin polarizations are 
FM(AF), there are two kinds of collinear configura- 
tions(see Fig.Efa)-(d)), which are denoted by AF-FM or 
FM-FM(AF-AF or FM-AF), according to whether the 
coupling between the left and right edges on the same 
layers is AF or FM. Particularly, when both edges are 
spin AF polarized, the interedge coupling is expected 
to be rather small for a relatively wide ribbon, leading 
to the nearly degenerate AF-AF and FM-AF states 5 ^. 
There also exists the possibility that the edges has differ- 
ent spin polarizations. The state denoted by figure Elm) 
is called " Mixed" , since it is spin FM polarized at the left 
edge but AF polarized at the right one. These five states 
fall into the first category of the magnetic states. Upon 
further doping at finite bias, the right edge of the up- 
per layer firstly becomes depolarized with the other spin 
polarizations qualitatively unchanged. These states are 
labeled by us by following the denotations of the states of 
the first category a symbol "-II". Correspondingly four 
states(shown in figure[2je)-(h)) are found to be stabilized, 
and are classified as the second category. Further dop- 
ing and increasing V will further depolarize the system. 
The third category consists of all the possible states with 
the upper layer fully spin depolarized. Since in this cate- 
gory, only the lower layer are spin polarized, analogously 
to MZGNR case, one can denote simply the magnetic 
states as AF or FM(see figure Eli), (j)), according to the 
nature of their interedge coupling. The fourth category 
includes two states where only one edge of the lower layer 
is spin polarized, as shown in figure EJk)-(l). These two 
states are denoted by FR and FL respectively. The fifth 
category consists of only one state(see figureEtn)), which 
lies near V = at finite doping and is denoted by FM- 
II, possessing the FM coupled edge polarizations located 
at the left edge of the upper layer and right edge of the 
lower layer. 

At half-filling, these exist two magnetic phases where 




(m) Mixed 



(n) FM-II 



FIG. 2: (Color online) Fourteen types of collinear magnetic 
structures for a BZGNR which has a width W — 10 and is pe- 
riodic along the vertical direction. The black(green) hexagons 
denote the upper(lower) layer lattices while the red(blue) ar- 
rows represent the spin-up(spin-down) polarizations at the 
lattice sites one the edges. 



the AF-FM phase is stabilized at lower V while the FM- 
AF one is stabilized instead at larger V. When V is larger 
than the critical value V c rs 0.105£, the BZGNR system 
will undergo a phase transition from the FM-AF state 
to a paramagnetic(PM) state. As mentioned before, the 
model Hamiltonian at half-filling is invariant under PI. 
It is found that all the above three states at half-filling 
preserve the PI symmetry and so each of them has the 
following symmetry properties: 
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are the electron densities and 




S 0.4 1*™'" 



0.00 



0.04 AF-FM 0.08 0.12 

V/t 



0.16 



0.20 



FIG. 3: (Color online) The magnetic phase diagram of a 
BZGNR with width W — 10: doping value Sn versus the 
bias voltage V/t, where the effective doping value Sn = 
(1 — n) x 2W is measured as the density departure from half- 
filling per unit cell of the bilayer ribbon, with n the average 
electron number per site. The dashed line(olive) marks the 
states in which the lower layer is at half-filling. The HM 
regimes are denoted and distinguished by the grid, cross- 
ing and vertical dashed lines, which represent that the cor- 
responding HM states possess distinct values of the quan- 
tized electric conductance 2e 2 /h, 3e 2 /h and 4e 2 //i respec- 
tively, along a perfect ribbon. 



edge spin polarizations at the sites on the four edges 
of the two layers. This PI symmetry is even preserved 
near the transition point between the AF-FM and FM- 
AF states, while it is completely lost when the system is 
away from half- filling (See figure 0}. 

On the other hand, in the absence of a gate bias, i.e., at 
V = 0, the model system is inversion symmetric, which 
means that the ground state might be invariant under 
inversion /. We find that upon increasing Sn, the inver- 
sion symmetric and inversion asymmetric states occur 
alternatively. Concretely, the ground state is found to 
be inversion symmetric for all the FM-FM states(0.05 < 
Sn < 0.325) and all the FM-II states(0.38 < Sn < 0.475) 
at V = 0, while found to be inversion asymmetric other- 
wise. There exist a regime of symmetric NC canted states 
located at < Sn < 0.02, which is <r v I symmetric apart 
from a global spin rotation, with a v the mirror reflection 
with respect to the plane normal to the zigzag direction. 
All the above inversion symmetric collinear states satisfy: 
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FIG. 4: (Color online)Amplitudes of the four edge polariza- 
tions as functions of the bias V/t (a) at half-filling, and (b) at 
Sn — 0.25. Superscripts and subscripts in m^jjL denote edge 
polarizations of the edge sites of the up(down) layer and at 
the left(right) edge, respectively. 



where the relative angle between m^ and m^ varies from 
to 33°. 

Once V > V c , sufficient electrons will be transferred 
from the upper to the lower layer, making the upper layer 
a PM state without polarizations at both edges, regard- 
less of the doping level. Since the upper layer is fully 
depolarized, the magnetism of the BZGNR is dominated 
by the lower layer and then the system can be viewed 
as one possessing approximate particle-hole symmetry at 
half-filling of the lower layer. As a consequence, the mag- 
netic states together with the PM ones are found to dis- 
tribute nearly symmetrically about the half-filling line of 
the low layer (the dashed line in figure [3]) in the phase 
diagram. The other consequence is that starting from an 
AF state, upon doping, the system undergoes transitions 
one by one to the FM, to the FR, and to the PM states 
respectively, analogous to that in MZGN R 21 i 45 . Further- 
more, except for the NC canted states, all phases located 
at this dashed line, including the AF, the AF-FM-II and 
the Mixed states, have AF couplings between the left and 
right edges on the lower layer, which is consistent with 
that in MZGNR^Sa^. 

The topological structure of the phase diagram is 
nearly unchanged for different ribbon width W as long as 
it is less than 20. For a wider ribbon, we note that since 
the interedge coupling is more negligible, especially for 
the case where one of the edges is spin AF polarized—, 
the FM-AF and AF-AF states (or the FM-AF-II and AF- 
AF-II states) become competing and in certain regime of 
phase diagram, they are nearly degenerate. 



IV. THE CANTED STATES AND THE HM 

STATES 



m, 



m, 



m, 



m 



n- 



while the symmetric coplanar NC states satisfy: 
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Now let's focus on the magnetic properties of the NC 
W canted states. It can be seen from the phase diagram that 
the NC canted states exist in the low-doping and small- 
bias regime, which is much larger, compared with that 
in MZGN R 21 ! 45 . So the magnetic interlayer interactions 
really alter the character of the magnetic states. All the 
(5) NC canted states are found to be coplanar and so can 
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FIG. 5: (Color online)The left panels (a), (c) and (e) ex- 
hibit the relative orientation angles between the spin polariza- 
tion directions, while the right panels (b), (d) and (f) exhibit 
the magnitudes of the four edge polarizations as functions of 
doping, at different biased voltages V . Two representative 
band structures of the NC canted states are shown in (g) 
at V/t — 0, 6n — 0.015 for a symmetric one, and in (h) at 
V/t = 0.04, Sn = 0.175 for an asymmetric one. 8 u (6 d ) are 
the relative angles between the left and right edge spins on 
the upper(lower) layer, and 6l(8r) are the relative angles be- 
tween the left(right) edge spins on the upper ad lower layers. 
l m ifiol are defined the same as that in figure [4] 



be characterized by the magnitudes of and the relative 
orientation angles between edge polarizations. We denote 
6 u (8 d ) as the relative angles between the left and right 
edge spins on the upper(lower) layer, and 9^(9 R ) as the 
relative angles between the left(right) edge spins on the 
upper ad lower layers. 

These NC canted states are always found to be metal- 
lic, as can be seen by the band fillings of the two typical 
states shown in figure EIg)and (h). It can be seen that 
among the eight split edge bands due to the bias voltage 
as well as the Coulomb interactions, one(two) band(s) 
is(are) partly filled for the first(second) case. Therefore, 
dependent on doping, the spectrum is characterized by 
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FIG. 6: (Color online) Band structures for the seven typical 
HM states for the BZGNR with width W = 10: (a) V/t = 
0.06 and Sn = 0.015; (b) V/t = 0.09 and Sn = 0.025; (c) 
V/t = 0.09 and Sn = 0.125; (d) V/t = 0.09 and Sn = 0.1625; 
(e) V/t = 0.09 and Sn = 0.2625; (f) V/t = 0.14 and Sn = 
0.1375; (g) V/t = 0.14 and Sn = 0.1125. 



four or eight counterpropagating current-carrying states 
at the Fermi level, leading to a quantized electric con- 
ductance 2e 2 /h or 4e 2 /h along the ribbon with perfect 
edge, respectively. 

Three representative doping processes related to the 
NC canted states are shown in figure E£a),(c),(e). For 
the unbiased BZGNR 7 = 0, the BZGNR at half-filling is 
the AF-FM insulating stated. Upon doping the ground 
state undergoes a continuous phase transition to a a v I 
symmetric NC canted state, which is characterized by the 
relations 9 U = 9 d , L = 6 R . When 0.02 < Sn < 0.052, 
the NC ground states becomes asymmetric, i.e., 9 U ^ 9 d 
and 9 L ^ 9 R . 

figure [5jb),(d),(f) give the edge magnetizations of a 
BZGNR as functions of doping for different V, which 
clearly exhibit the depolarization process upon doping. 
In particular, the latter two show the successive depolar- 



ization processes via the first-order transitions for cases 
at finite bias. The right edge of the upper layer becomes 
firstly depolarized, then the left edge of the upper layer, 
the left edge of the lower layer, and finally the right edge 
of the lower layer. This is interpreted as follows. When a 
BZGNR is doped away from half-filling, a finite gate bias 
will transfer electrons from the upper to the lower layer 
and then drive the upper layer further away from half- 
filling while drive the lower layer closer to half-filling, 
making this layer easier to be spin-polarized. On the 
other hand, the left(right) edge sites of the upper(lower) 
layer are most easily spin polarized since they are actu- 
ally the leftmost (rightmost) sites of the BZGNR, which 
are closer to magnetic instability. 

In analogy with MZGNR, the half-metallicity can be 
expected in BZGNR. The regimes where the HM states 
exist are clearly demonstrated in the phase diagram. Dif- 
ferent HM states can be distinguished by their magnetic 
orders and the values of their quantized electrical conduc- 
tances along the ribbon, the latter of which can be de- 
termined by counting the number of the current-carrying 
states at the Fermi level, i.e., the number of the crossing 
points between the band spectrum and the chemical po- 
tential. From the phase diagram, we see that all the HM 
states form a connected regime below the half-filling line 
of the lower layer, indicating that the HM states only ex- 
ist where the doping level of the lower layer is more than 
half-filling. We find seven kinds of the HM states in five 
magnetic phases(see figure [5]): The AF-FM phase at fi- 
nite doping is HM with a quantized electric conductance 
2e 2 //i; The FM-AF phase is HM and has a quantized 
electric conductance 2e 2 /h at low doping but 4e 2 //i at 
larger doping. The FM-AF-II phase are HM one charac- 
terized by a quantized electric conductance 3e 2 //i. Two 
parts of the FM phase are HM with quantized electric 
conductance 2e 2 /h and 4e 2 //i respectively. Finally, one 
part of the FR phase is HM, possessing a quantized elec- 
tric conductance 3e 2 //i. Thus for a BZGNR at a fixed 



doping, one can achieve phase transitions between dif- 
ferent HM states by tuning the bias voltage, which may 
have great applications in graphene manipulation. 



V. SUMMARY 

In summary, we have investigated the possibility of the 
noncollinear spin canted states and half-metallicity in a 
bilayer zigzag nanoribbon under a bias voltage. Due to 
the interlayer and intralayer couplings between the edge 
polarizations, the magnetic states are classified and the 
phase diagram is given. At a small but finite bias voltage, 
upon doping away from half-filling, the bilayer ribbon will 
be spin depolarized step by step: firstly the right-edge of 
the upper layer, then the left-edge of the upper layer, the 
left-edge of the lower layer, and finally the right-edge of 
the lower layer. Both symmetric and asymmetric non- 
collinear spin canted states are found to be stabilized 
as metallic ones in the low-doping and small-bias regime. 
Seven kinds of the HM states are found, each of which has 
distinct magnetic structure or different quantized elec- 
trical conductance along the perfect ribbon. Since the 
HM properties are important in spin-related transport, 
these results are expected to have potential applications 
in bilayer-graphene-based spintronics. 
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